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We study the pure SU(3) gauge theory in 2+1 dimensions on the lattice using 't Hooft's twisted 
boundary conditions to force non-vanishing center flux through the finite volume. In this way we 
measure the free energy of spacelike center vortices as an order parameter for the deconfinement 
transition. The transition is of 2"'' order in the universahty class of the 2d 3-state Potts model, 
which is self-dual. This self-duality can be observed directly in the SU(3) gauge theory, and it 
can be exploited to extract critical coupUngs with high precision in rather small volumes. We 
furthermore obtain estimates for critical exponents and the critical temperature in units of the 
dimensionful continuum coupling. Finally, we also apply our methods to the (2+l)ii SU(4) gauge 
theory which was previously found to have a weak order transition. We nevertheless observe 
at least approximate q = A Potts scaling at length scales corresponding to the lattice sizes used in 
our simulations. 
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1. Introduction 

The motivation of our study is to see how much we can learn about the deconfinement tran- 
sitions in pure SU(A^) gauge theories at finite temperature from universality and scaling. As these 
concepts require second order phase transitions, in 3+1 dimensions we are only left with SU(2), 
where spatial center vortex sheets share their universal behavior with interfaces in the 3d Ising 
model. This has been studied in [1]. In 2+1 dimensions on the other hand, we can apply analogous 
methods to both SU(2) and SU(3) which then exhibit a 2"'' order deconfinement transition. The 
latter is in the universality class of the 2d Ising model for which by far the largest pool of exact 
results is available. This was exploited in [2] for a high precision determination of critical cou- 
plings and temperature, and in [3] for accurate finite-size scaling, a reflection of self-duality and 
a precision determination of the behavior of the continuum sting tension and its dual around the 
phase transition. Here, we report on first results from applying the same methods to SU(3) and 
gradually also to SU(4) in 2+1 dimensions. 

For SU(3), the dimensionally reduced spin model with the same global symmetry and the uni- 
versal properties of the Z3 center-symmetry breaking deconfinement transition is the 3-state Potts 
model. For SU(4), the Z4 center symmetry alone does not uniquely specify the effective spin model 
to describe the dynamics of Polyakov loops. SU(4) is a rank-thi-ee group and has three fundamen- 
tal representations, 4, 4 and 6. So even the simplest effective Polyakov-loop model will consist 
of two distinct real terms, with nearest neighbor couplings between loops in 4/4 representations 
and between loops in the 6 representation [4]. Depending on the relative weight between the two, 
the corresponding spin model could be any of the Z4-symmetric Ashkin-Teller models with three 
energy levels per link and continuously varying critical exponents between the q = 4 Potts model 
class and that of the planar or vector Potts model which corresponds to two non-interacting Ising 
models in this case. Even though the more recent studies of the {2+l)d SU(4) gauge theory [5, 6, 7] 
indicate that the transition is weakly P' order, we do find, at least approximately at the length scales 
corresponding to our spatial lattice volumes, a universal scaling which seems closest to the standard 
q = 4 Potts case. Because this might not seem very natural, it would be interesting to explain that. 

One special feature of the ^-state Potts models in 2 dimensions is that they are self-dual for all 
q, no matter whether they have {q> 4) or 2"'^ (q < 4) order transitions. With a 2"'' order transition 
and scaling, this self-duality is reflected in the (2+l)<i gauge theory: the spatial center-vortex free 
energies are miiTor images around criticality of those of the confining electric fluxes [3]. Here we 
verify this explicitly for SU(3), and show how it can be used to remove the leading (universal) 
finite-size corrections in the determination of critical couplings from finite volume extrapolations. 

2. Concepts and Methods 

't Hooft's twisted boundary conditions, center vortices and electric fluxes: In a theory with- 
out matter fields where the gauge fields represent the center of the gauge group trivially, the bound- 
ary conditions in a finite Euclidean l/T x L'' volume are only fixed up to center elements giving 
rise to N'^ gauge-inequivalent boundary conditions for a pure SU(A'^) gauge theory in d+\ dimen- 
sions. These twisted boundary conditions [8] can be classified either as magnetic twists defined in 
purely spatial planes or as temporal twists in the planes oriented along the Euclidean time direction. 
The latter are labeled by a vector k G Zfj. In the following we will only consider temporal twists 
because the magnetic twists are inelevant for the deconfinement transition. 
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Temporal twist introduces spatial center vortices whose free energies provide order parame- 
ters for the deconfinement transition. These vortex free energies (per temperature T) are de- 
fined as ratios Rk(k) = Zk(k) /Zk(v) = of partition functions Zk(k) with temporal twist k 
over the periodic ensemble Z)t(O). Analogously, one defines the electric flux free energies via 
Re{e) = Ze(e) jZgi^) = e^^'^'^\ These describe gauge-invariant color-singlet free energies of static 
fundamental charges at some point x with mirror anti-charges in a neighboring volume a.tx + Le 
along the direction of the flux e relative to the no-flux ensemble Ze(0) = J^^Zii(k), which is an 
enlarged ensemble corresponding to fluctuating temporal twists, see [1, 3]. The electric flux and 
spatial center vortex partition functions are related by a t/-dimensional Z^ Fourier transform, 

R,ie) = i(tr(P(x)pt(x + eL))>„„.,^^ = {Hie^^^''/^ R,(k)) /^^^R,{k) . (2.1) 

Universality and self-duality: By the Svetitsky-Yaffe conjecture, a. d + I dimensional gauge 
theory with second order deconfinement transition has the same universal properties as a J dimen- 
sional spin model with the same global Zf^ symmetry [9]. The Polyakov loop correlators of the 
gauge theory near criticality behave in the same way as those of spins in the spin model. Spatial 
center vortices con^espond to spin interfaces, which are frustrations where the coupling of adjacent 
spins favors cyclically shifted spin states rather than parallel ones for the usual feiTomagnetic cou- 
plings. Consequently, the center-vortex free energies F^ show the universal behavior of interface 
free energies. For the 2-1-1 dimensional SU(3) gauge theory the corresponding spin model is the 2d 
3-state Potts model. Like all <7-state Potts models in 2 dimensions, it is self-dual. 

Kramers-Wannier duality is of course a very well-known concept from statistical physics [11]. 
It provides exact maps between the spin systems and their dual theories in terms of disorder vari- 
ables on the dual lattice. In 2 dimensions, just as for the Ising model (q = 2) these dual theories 
of the <7-state Potts models ai^e again <7-state Potts models, but at a dual temperature T which is 
swapped around criticality at Tc as compared to the original model. Duality transformations in a 
finite volume do not preserve boundary conditions, however. Periodic boundary conditions on one 
side generally correspond to fluctuating boundary conditions on the other [12]. This was explicitly 
demonstrated for the duality between the 3d Ising and the Z2-gauge model in [13]. The exact finite- 
volume duality transformation for the 2d <7-state Potts models is given in [3]. Here it suffices to 
note that its structure is precisely that in (2.1) with q = N. It expresses the partition function of the 
dual ^-state Potts model with certain set boundary conditions at a temperature T as a 2t/ Z^, -Fourier 
transform over Potts models with all possibilities of cyclically shifted boundary conditions at tem- 
perature T. In the {2+l)d gauge theory, the temperature is the same on both sides of the ZAr-Fourier 
transform (2.1). But within the universal scaling window around a 2"'' order phase transition, as a 
consequence of the self-duality of the spin model, the free energies of spatial center vortices and 
those of the confining electric fluxes are mirror images of one another around Tc. 

Numerical procedure: To implement 't Hooft's twists on the lattice we multiply a stack of pla- 
quettes by the corresponding center element z&Z^ so as to fix the con^esponding amount of center 
flux through the planes with twisted boundary conditions. For the temporal ^-twists this introduces 
thin spatial center vortices perpendicular to k which separate regions where fundamental Polyakov 
loops differ by the center phase z = e^^''^^^. They are thus like spin interfaces in the Potts model. 
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Figure 1: Self-duality in SU(3) (A^, = 24, Nt = 2) 



We use the snake algorithm [14] to compute the ratios Rk(k) = Zi^(k) /Zj^{0) via products of 
ratios of partition functions which differ by only one phase rotated plaquette to successively form a 
vortex that wraps around the spatial volume perpendicular to k. We use the standard Wilson gauge 
action and a heatbath algorithm supplemented by a number of overrelaxation steps. SUCA'^) updates 
are implemented as updates of SU(2) subgroups in the usual way [15]. 

3. Results 

Self-duality in SU(3): As for SU{2) in 2+1 dimensions [3] the self-duality of the 3-state Potts 
model is reflected in SU(3). The spatial center vortex and electric flux partition functions are related 
by a Z3-Fourier transform as in (2. 1) whose structure is the same as that of the finite volume duality 
transformation of the Potts model. Its self-duality therefore implies that electric fluxes on one side 
of the phase transition should resemble center vortex ensembles on the other. 

Swapping the temperature T -f-)- T in the spin model near criticality amounts tox ^ —x where 
X = ±L/^± is the finite-size scaling variable given by the ratio of the finite size L over the correla- 
tion lengths ^± = f±{±t)^^, with for x > and t >0 above Tc, and for x,t <0, below. For 
the gauge theory we use x = ±TcL (ib?)^ocL/(^* and control the reduced temperature t = T/Tc—l 
by changing the lattice coupling. Within the universal scaling window, as functions of x, we should 
then find that Rk{x) = Re{—x) for matching pairs of twist k and flux e. This is indeed the case also 
for SU(3), as demonstrated in Fig. 1, where we compare the ratios of partition functions and R^ 
for one unit of temporal twist and one unit of electric flux, i.e., both k and e either (1,0) or (0, 1), 
over the phase transition. 

Critical couplings: There is a long history of methods to extract critical couplings or temper- 
atures from simulations in finite volumes, going back to using pairwise intersections of Binder 
cummulants on sucessively larger lattices [16]. Hasenbusch later demonstrated that the ratios of 
partition functions with different boundary conditions could be used in the same way to obtain a 
much more rapid convergence with very good estimates already from rather small lattices [17]. At 
criticality, these ratios tend to universal values < /?c < 1 in the thermodynamic limit. In [2] it was 
therefore shown how to obtain critical couplings for gauge theories from intersecting the ratios R^ 
of finite volume partition functions with these universal fixed points, once their values are known. 
For {2+l)d SU(2) this led to an even faster convergence than their pairwise intersections. For the 
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Figure 2: Convergence to j3,. (Nt = 4) fi-omi?i(j3) =i?e(j3) and7?^(j3) in SU(2) (left) and SU(3) (right). 

2d Potts models with 2"'* order transition, i.e. , for ^ = 2, 3 and 4, the universal numbers have 
been obtained exactly, in terms of Jacobi theta functions, for all cyclic boundary conditions (with 
m, « = 0, 1 , . . . g — 1) in [10]. For ^ = 3 on a symmetric lattice they are, 

4' °) =/?P'°' =0.30499982... , and 4''^' =/?P^ =/?^^' = 0.19500018 .. . . (3.1) 

Using a finite-size scaling ansatz for the vortex ensemble ratios 7?^- around criticality of the form 

Rk{l5)=Rc + b{l5- j8,) nI^" + cN-'" + ---, (3.2) 
we define pseudo-critical couplings l5c{Ns,Nt) in a finite volume by requiring that Rk{P) = Rc, 

MK,N,) = MN,) - {c/b)N;^''^"''^ + • • • . (3.3) 

These extrapolate to ^c{Nt) from large spatial lattice sizes A'^v at fixed numbers of time slices A',. As 
a byproduct this method gives numerical estimates of the correction to scaling exponent ft). 

With self-duality, however, there is a 
yet more efficient method to determine jSc 
[3]. This is based on the simple observation 
that one must then have Re{^) = Rk{^) for 
like e and ^ at /3 = /3c. In fact, one can eas- 
ily convince oneself that with self-duality, 

R,{l5)=R,-b{l5- ^,)nI" +cN;"' + • • • , 

with the same coefficients b and c as in (3.2) 

Therefore, the leading finite-size corrections to jSc when defined by Re = Rk cancel. At critical- 
ity, ReiPc) = Rkil^c) = Rc + cN^"' H , SO the leading corrections only move the intersection point 

vertically without shifting the so defined critical coupling. The gain is illustrated for both SU(2) 
and SU(3) in Fig. 3. When intersecting R^ = R^, we form weighted means from sufficiently large 
aspect ratios Ns/Nt where we assume that the estimates have converged within errors. These ai^e 
compared to the extrapolated values from intersecting Rk{P) with the universal value Rc in Tab. 1. 
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8.15297(57) 
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Table 1: SU(3) critical couplings from self-duality 
(weighted means), intersection with the universal value 
(extrapolated), and literature values from ^^[7], *[18]. 
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Critical temperature and correlation length exponent v: In 2+1 dimensions the critical cou- 
plings grow linearly with Nt to leading order at large A^,. The slope is given by the critical temper- 
ature in units of the dimensionful continuum coupling, ^c{Nt) / {2Nc) = {Tc/gl)Nt H . 

From our values for Nt = 4, 6 and 8 we then obtain Tc/g^ = 0.5475(3) corresponding to 
Tel \/o = 0.9938(9) with a zero temperature string tension ^/o /g\ = 0.5509(4) from a weighted 
average of the four values in [19]. This is consistent with Td y/o = 0.9994(40) from [7]. 

Moreover, because the spatial center vortex free energies F^. for sufficiently large L depend 

only on V-/^t, and t oc (j8 - j3c), when expanding Fk{fi) = - ln/?c + <i(A^v)(jS - jS,) H , we can 

expect the slope at j8c to behave as d{Ns) ~ A^,^''^- The result from fitting our slopes for SU(3) with 
Nt=4 then gives v = 0.82(4) as compared to v = 5/6 f« 0.833 for the 2d 3-state Potts model. 

Results for SU(4): The present conclusion from a sequence of studies of the {2+\)d SU(4) gauge 
theory [5, 7, 6] is that the deconfinement transition is weakly order. Especially the detailed 
analysis in [6] was consistent with first order volume scaling laws. Here we assess to what extent 
Potts model scaling describes the transition, at least approximately, and whether we find indications 
of where our methods start to fail as we go to larger and larger volumes. 
First, we extract critical couplings from i 



the pairwise intersections of the Fyt's for pairs 
of lattices with Ns ratios of 2: 1 . This method 
is independent of Potts scaling and yields 
j3c = 26.283(9) for A^, = 4. Then we com- 
pare the so extrapolated value to that ob- 
tained from intersecting the F/t's with the 4- 
state Potts universal value from [10]. The 
latter has a smaller enor because we have 
more points to fit; it gives the consistent value 
j3c = 26.294(2). Pseudo-critical couplings 
and fits for each method are shown in Fig. 3. 
Both extrapolated values are consistent with 

j8c = 26.228(75) from [7] but deviate with some significance from jS^ = 26.251(16) given in [6], 
where first order scahng was assumed in the infinite volume extrapolation of the critical coupling. 
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Figure 4: Check of Potts model scaling (A^, = 4) 
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If we furthermore extract a critical exponent v from the slopes d{Ns) of Fk{P ) at jS^ as before, 
see Fig. 3, we obtain v = 0.60(2) from our data for Nt = 4 with j8c = 26.283. Some systematic 
uncertainty arises from what precise value is used here, however. Repeating the analysis for differ- 
ent values in the one-a interval around jSc = 26.283, we find V = 0.59(5). Of all the Z4-symmetric 
Ashkin-Teller models, which have continuous v S [2/3, 1], this seems to be at best consistent with 
the lower bound v = 2/3 for the q = 4 Potts model, in agreement with the earlier conclusion in 
[5]. It is a general trend of our method, observed also for SU(2) and SU(3), that it underestimates 
the critical exponent due to subleading finite-size effects, however. Finally, our present A/, = 4 data 
with spatial lattice sizes up to A'^, = 80 shows reasonably good Potts scaling as seen in Fig. 4 where 
we plot the center-vortex free energy Fk{x) over the scaling variable x = ±TcL{±ty with v = 2/3. 

4. Conclusions 

We have studied the deconfinement transition in the pure SU(3) gauge theory in 2-1-1 dimen- 
sions on the lattice. Using 't Hooft's twisted boundary conditions we have measured center-vortex 
free energies and demonstrated that the self-duality of the associated Potts model is directly re- 
flected in SU(3): the free energies of the confining electric fluxes ai^e mirror images ai^ound Tc of 
those of spatial center vortices. We demonstrated how this can be exploited to remove the leading 
finite-size connections in the determination of critical couplings from numerical simulations. We 
do not yet have the data necessary to compute electric fluxes and to test self-duality analogously in 
SU(4), but our available data does not show any significant violations of the q = 4 Potts scaling. 

Acknowledgements: This work was supported by the Helmholtz International Center for FAIR 
within the LOEWE program of the State of Hesse, the Helmholtz Association Grant VH-NG-332, 
and the European Commission, FP7-PEOPLE-2009-RG No. 249203. Simulations were performed 
on the high-performance computing facilities of eResearch SA, South Australia. 
References 

[1] Ph. de Forcrand and L. von Smekal, Phys. Rev. D 66 (2002) 01 1504; Nucl. Phys. (PS) 106 (2002) 619. 

[2] S. Edwards and L. von Smekal, Phys. Lett. B 681 (2009) 484. 

[3] L. von Smekal, S. R. Edwai'ds and N. Strodthoff PoS Lattice 2010 (2010) 292; 
N. Strodthoff, S. R. Edwards and L. von Smekal, in preparation. 

[4] C. Wozar et al, Phys. Rev. D 74 (2006) 1 14501; A. Wipf, private communication. 

[5] Ph. de Forcrand and O. Jahn, Nucl. Phys. Proc. Suppl. 129 (2004) 709. 

[6] K. Holland, M. Pepe and U. J. Wiese, JHEP 0802 (2008) 041. 

[7] J. Liddle and M. Teper, arXiv:0803.2128. 

[8] G. 't Hooft, Nucl. Phys. B 153 (1979) 141. 

[9] B. Svetitsky and L. G. Yaffe, Nucl. Phys. B 210 (1982) 423. 
[10] H. Park and M. den Nijs, Phys. Rev. B 38 (1988) 565. 
[11] R. Savit, Rev. Mod. Phys. 52 (1980) 453. 

[12] C. Gruber, A. Hintermann and D. Merlini, Group Analysis of Classical Lattice Systems, Springer, 

Berlin - Heidelberg 1977; M. Caselle and M. Hasenbusch, private communications. 
[13] M. Caselle, M. Hasenbusch, P Provero and K. Zarembo, Nucl. Phys. B 623 (2002) 474. 
[14] Ph. de Forcrand, M. D'Eha and M. Pepe, Phys. Rev. Lett. 86 (2001) 1438. 
[15] N. Cabibbo and E. Marinari, Phys. Lett. B 119 (1982) 387. 

[16] K. Binder, Z. Phys. B 43 (1981) 119; K. Binder and E. Luijten, Phys. Rept. 344 (2001) 179. 

[17] M. Hasenbusch, Physica A 197 (1993) 423. 

[18] J. Engels, et al., Nucl. Phys. Proc. Suppl. 53 (1997) 420. 

[19] B. Bringoltz and M. Teper, Phys. Lett. B 645 (2007) 383. 



7 



